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THE  INFLUENCE  OF  TIP  MASS  OFFSET  ON 
THE  STABILITY  OF  BECK'S  COLUMN 

In  this  paper,  the  stability  of  a slender  cantilever  carrying  a tip 
mass  at  its  free  end  and  subjected  there  to  a follower  force  is  investigated. 
The  centroid  of  the  tip  mass  is  offset  from  the  free  end  of  the  beam  and  is 
located  along  its  extended  axis.  The  associated  boundary  value  problem  is 
solved  and  the  exact  frequency  equation  is  derived.  The  frequency  equation 
is  solved  numerically  for  the  case  in  which  both  the  beam  and  the  tip  mass 
have  circular  cross  sections.  The  numerical  computations  indicate  that  the 
system  loses  stability  only  through  flutter.  The  variation  of  the  values 
; of  the  critical  flutter  load  Qcr  with  the  tip  mass  offset  parameter  £ is 

shown  graphically  for  four  values  of  the  tip  mass  density  to  beam  density 
ratio  p.  These  calculations  reveal  that,  at  sufficiently  small  values  of 
£,  Qcr  decreases  sharply  for  increasing  values  of  p.  For  values  of  £ 
sufficiently  large,  however,  the  situation  is  reversed  as  the  value  of  Qcr 
increases  with  increasing  p. 

1 . INTRODUCTION 

An  elastic  cantilever  of  length  SI  and  density  p that  is  subjected  to  a 
compressive  follower  force  of  magnitude  P applied  at  its  free  end  is  known 
as  Beck's  column  [1].  Pfluger  [2]  has  investigated  the  influence  of  the 
transverse  inertia  of  a tip  mass  on  the  value  of  the  critical  flutter  load 
of  Beck's  column,  and  Walter  and  Levinson  [3]  and  Anderson  [4]  have  con- 
sidered closely  related  problems  with  the  inclusion  of  the  influence  of 
the  rotatory  inertia  of  the  tip  mass.  If  the  rotatory  inertia  of  the  tip 
mass  is  negligible,  its  transverse  inertia  generally  tends  to  reduce  the 
. value  of  the  critical  flutter  load.  However,  if  the  mass  of  the  column  is 
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sufficiently  small  relative  to  that  of  the  body  attached  at  its  tip,  the 

value  of  the  critical  flutter  load  can  be  slightly  greater  than  Beck's 

value  of  P = 20.05  EI/£2. 
cr 

In  all  the  studies  mentioned  above,  it  was  assumed  that  the  mass  of  the 
attached  body  was  concentrated  at  the  free  end  of  the  column.  As  a further 
step  in  studying  the  influence  of  a tip  mass  on  the  stability  of  Beck's 
column,  one  may  consider  the  system  in  which  the  centroid  of  the  tip  mass  is 
offset  from  the  point  of  attachment  along  the  extended  axis  of  the  column. 

Only  recently  the  natural  frequencies  of  free  vibration  of  a uniform  unloaded 
cantilever  carrying  such  a tip  mass  have  been  reported  by  Bhat  and  Wagner  [5], 
Bhat  and  Kulkarni  [6] , and  Flax  [7] . 

2.  STATEMENT  OF  THE  PROBLEM 

The  differential  equation  of  motion  and  the  boundary  conditions  for  the 
system  consisting  of  Beck's  column  carrying  at  its  free  end  a tip  mass  whose 
centroid  is  offset  from  the  point  of  attachment  (see  Figure  1)  are 

HI  + p l!”  + PA  — = 0,  0 < x.  < Z , (1) 

Bxj1*  3x^2  3t2 


dw 

W(0,t)  = — (0,t)  = 0 , 
3xi 


(2) 


3zw 


33w 


El  — (£,t)  = - (JQ  + me2)  (Z,t)  - mc^-  (Z,t)  , 


3x^2 


3x^3t 


32_w 
3t2 


(3) 


El  — (£,t)  = m — (Z,t)  * me  (Z, t)  , 

3xj3  3t2  3xj3t2 


(4) 


where  w(xj,t)  denotes  the  transverse  deflection  of  the  column,  El  its 
flexural  rigidity,  P is  the  magnitude  of  the  applied  force,  p the  density 
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of  the  beam,  A its  cross-sectional  area,  l its  length,  JQ  the  moment  of 
inertia  of  the  tip  mass,  and  c the  distance  from  the  end  of  the  beam  to  the 
centroid  of  the  tip  mass.  If  one  sets  c = 0 in  equations  (3)  and  (4),  then 
the  given  boundary  value  problem  reduces  to  the  problem  that  is  the  adjoint 
of  the  one  solved  in  reference  [3],  whereas,  if  both  c and  JQ  are  equated 
to  zero,  Pfluger's  problem  [2]  is  obtained. 

It  is  expedient  to  express  equations  (1)  to  (4)  in  dimensionless  form 
by  making  the  changes  of  independent  variables  Xj  = 2.x,  0 < x < 1,  t = gi, 
and  by  introducing  the  following  parameters: 


g>  = p^:  Q = n\  Y 

EI  El 


_2_  u = JL_ 
pA23  pA2 


a = — . 


As  a consequence  of  these  changes,  equations  (1)  to  (4)  can  now  be  expressed 
as  follows: 

w""  + Qw"  + w = 0,  0 < x < 1,  0 < t,  (6) 

W(0,T)  = w'(O.T)  = 0,  (7) 

w"(1,t)  + (y  + ya2)w'(l,T)  + uaw(l,x)  = 0,  (8) 

w"» (l,r)  - uw(l,T)  - paw'(l,T)  = 0,  (9) 

where  primes  and  dots  denote  derivatives  with  respect  to  x and  t,  respectively. 

3.  THE  EXACT  SOLUTION 

The  solution  of  the  partial  differential  equation  (6)  is  next  assumed  in 
the  form 


w(x,x)  = y(x)e  , (1 

1/2 

where  i = (-1)  and  a denotes  the  dimensionless  natural  frequency  of 
vibration  of  the  beam.  Substitution  of  equation  (10)  into  equations  (6)  to 
(9)  yields  the  following  non-self-adjoint  eigenvalue  problem: 
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and 


y""(x)  + Qy"(x)  - w2y(x)  = 0, 


0 < x < 1, 


(ID 


(12) 

(13) 

(14) 

(15) 

(16) 

(12) 

(17) 

where 

= 1«  aj2  = a13  = a14  = 

a2i  = °»  a22  = Al*  a23  = °»  a24  = X2* 

a31  = ~(^i2  + w2ua)cosA1  + + ya^sinAj, 

a32  = ”^i2  + w2ua)sinXj  - A^fy  + ya2)cosXj, 

a33  = (^2*  ~ a)2lja)cosh^2  ' ^2a)2^Y  + lia2)sinh^2’ 

a34  = C^2  ” a,2yot)sin^^2  " ^2t°2^  + pa2)coshX2» 

a41  = AjfAj^  - w2ua)sinXj  + w2p  cosA^, 


y(0)  = y'(0)  = 0 

y"(l)  - a)2(Y  + ua2)y'(l)  - yoiw2y(l)  = 0, 
y"'  (1)  + yai2y(l)  + uaa)2y'(l)  = 0 . 

It  can  be  verified  that  the  solution  of  equation  (11)  has  the  form 
y(x)  = Aj  cosXjx  + A2  sinX^x  + Aj  coshX2x  + A4  sinhX2x  , 

where 

An  = (1//2)[(Q2  + 4oi2) 1/2  - C-DnQ] 1/2  , n = 1,2  . 
Substitution  of  equation  (15)  into  the  boundary  conditions  in  equations 
to  (14)  yields  the  following  system  of  homogeneous  algebraic  equations: 

4 

£ ajkAk  = °*  j = 1>2>3’4 
k=l 
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a42  = '^l^l2  " t»2ya)cosXi  + w2y  sinAj, 
a43  = ^2^2*  + w2ya)sinhA2  + w2y  coshX.,, 

a44  = ^2^2*  + u,ZlJa)c0S^^2  + w2P  san^^2  ’ 

The  system  of  equations  (17)  will  have  a non-trivial  solution  if  and 

only  if  the  determinant  of  the  coefficient  matrix  vanishes,  i.e., 

det(ajk)  = 0 . 

Expansion  of  this  determinant  yields  the  frequency  equation 

Q2  + 2w2  + 2uyw4  - a)  [Q(td2yY  - 1)  + ya(Q2  + 4w2)]sinA1  sinhX2  - 

- a)X2(X12  + X22)[y  + Aj2(y  + ya2)]sinX1  coshX2  + 

+ toAjfAj2  + X22)[y  - A22(y  + ya2)]cosXj  sinhX2  + 

+ 2w2(1  - yYu2)cosX1  coshX2  = 0.  (18) 

4 . NUMERICAL  RESULTS 

The  objective  of  the  numerical  computations  is  to  determine  the  depen- 
dence of  the  critical  load  parameter  Qcr  on  the  parameters  a,  y,  and  y.  which 
are  inter-dependent.  Hence,  it  is  necessary  to  examine  these  quantities 
somewhat  more  closely.  The  moment  of  inertia  parameter  Jg  is  defined  by 

J0  = / O42  + z22)dm,  dm  = P0dT,  (19) 

T 

where  t denotes  the  volume  of  the  tip  mass  and  its  density.  The  Zj^z^  - 

coordinate  axes  have  their  origin  at  the  centroid  of  the  tip  mass  as  shown 

in  Figure  2.  Since  dr  = dz^dz^Zj,  one  has,  under  the  assumption  that  the 
tip  mass  is  a prismatic  solid  whose  generators  are  parallel  to  the  Zj-axis, 


— 


i-  j— • x ... 


_ — 


! 


J0  = f pOc/  (cZ  + 3z22)dz2dz3, 


(20) 


where  Ag  denotes  the  cross-sectional  area  of  the  tip  mass.  To  proceed 
further,  one  must  assume  something  about  the  geometry  of  the  beam  and  of  the 
tip  mass.  For  the  sake  of  being  specific,  suppose  that  the  tip  mass  and 
beam  have  circular-cross  sections  of  radius  a and  b,  respectively,  where  it 
will  be  hypothesized  that  a > b.  Thus,  it  follows  from  equation  (20)  that 


J0  = l PQca  (3a2  + 4c2)  . 


(21) 


In  addition,  for  the  beam  of  circular  cross-section,  it  may  be  verified  that 


A = Trb2 


I = t b' 


(22) 


Therefore,  in  view  of  equations  (21),  (22),  and  m = 2pgivca2,  one  finds 
from  equation  (5) 


U = 


2pQa2c 

p£b2 


pQca2(3a2  + 4c2) 


6pb  £ 


2 o 3 


(23) 


Clearly,  for  a beam  of  given  dimensions  and  density,  the  values  of  y and  y 
can  be  changed  by  varying  a,  c,  and  pQ  either  together  or  independently. 
However,  because  the  essentially  new  parameter  here  is  c,  let  it  be  supposed 
that  Pq  and  a are  held  fixed  throughout  a given  sequence  of  computations 
while  c is  varied.  Thus,  the  objective  of  the  present  numerical  computations 
is  to  determine  Q as  a function  of  the  ratio  £ = 2c /£,  which  is  the  ratio 
of  the  length  of  the  tip  mass  to  the  length  of  the  beam.  Hence,  one  has, 
in  place  of  equation  (23) , 

P r 


u = — (~)2Z,  Y = r~  (£2  + 3s2) , 

2 p b 12 


(24) 


where  s = a/£. 
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Suppose  further  that,  again  for  the  purpose  of  being  specific,  b/Jl  = 1/20 
and  a/X,  = 1/10,  so  that  a/b  = 2.  Then  equation  (24)  becomes 

a = K/2,  y = 4pS,  Y = j p£(£2  + 3/100),  (25) 

where  p = pQ/p.  Upon  using  equation  (25),  equation  (18)  is  now  solved 
numerically  for  the  critical  load  Qcr  as  a function  of  the  tip  mass  length 
parameter  £ for  p = 1/10,  1/2,  1,  and  2.  The  numerical  procedure  consists 
of  selecting  a value  of  Q and  computing  the  corresponding  values  of  the 
frequencies  uj  for  the  first  two  modes  of  vibration.  The  value  of  Q is 
successively  increased  from  zero  until  the  first  and  second  frequencies 
coalesce  at  the  critical  value,  Qcr.  As  long  as  Q < Qcr»  these  frequencies 
are  real  numbers,  but  for  Q > Qcr  they  are  complex  conjugate  numbers.  In 
the  present  case,  loss  of  stability  through  divergence  is  not  possible.  The 
onset  of  flutter  is  signaled  at  Q = Qcr.  The  results  of  these  numerical 
computations  are  shown  in  Figure  3,  where  the  variation  of  Qcr/ir2  has  been 
plotted  versus  £ over  the  range  0 £ E,  £ 1 for  the  four  stated  values  of  p. 

This  figure  reveals  that,  for  the  values  of  p considered,  the  value  of 
Qcr  decreases  monotonically  with  increasing  E,  on  the  interval  0 < 5 £ 1. 
Moreover,  the  value  of  Qcr  initially  decreases  rapidly  from  Beck's  value 
of  Qcr  = 2.031ir2  as  E,  increases  from  zero.  Indeed,  at  least  for  sufficiently 
small  values  of  £,  the  rate  of  decrease  of  Qcr  increases  as  the  density 
ratio  p increases,  i.e.,  the  initial  slope  of  the  tangent  to  the  curve  in 
the  Qcr£-plane  becomes  increasingly  negative  as  the  value  of  £ is  increased. 

As  the  value  of  E,  is  further  increased,  the  slope  of  the  tangent  to  a given 
QcrC-curve  eventually  tends  to  increase  over  the  range  of  E,  considered. 
Therefore,  for  a given,  sufficiently  small  value  of  E,,  the  critical  load 
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decreases  as  p is  increased,  but  for  a given  sufficiently  large  \ 
the  critical  load  increases  as  p is  augmented. 
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Figure  3.  Variation  of  QCR//tt2  versus  £>  for  four 
values  of  p.  n 
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US  ARMY  MATERIALS  5 MECHANICS 

US  ARMY  MOBILITY  EQUIP  RSD  COMD 

RESEARCH  CENTER 

ATTN:  TECH  LIB 

1 

ATTN:,  TECH  LIB  - DRXMR-PL 

2 

FT  BELVOIR,  VA  22060 

WATERTOWN,  MASS  02172 

NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 
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EXTERNAL  DISTRIBUTION  LIST  (CONT) 


NO.  OF  NO.  OF 

COPIES  COPIES 


COMMANDER 

US  ARMY  RESEARCH  OFFICE 

P.O.  BOX  12211  1 

RESEARCH  TRIANGLE  PARK,  NC  27709 

COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB  1 

2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACT 


ATTN:  DRXPE-MT  1 

ROCK  ISLAND,  IL  61201 

CHIEF,  MATERIALS  BRANCH 

US  ARMY  R5S  GROUP,  EUR  1 

BOX  65,  FPO  N.Y.  09510 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 

ATTN:  CHIEF,  MAT  SCIENCE  DIV  1 

DAHLGREN,  VA  22448 

DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  DIR,  MECH  DIV  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 

COMMANDER 

WRIGHT- PATTERSON  AFB 

ATTN:  AFML/MXA  2 

OHIO  45433 

NASA  SCIENTIFIC  f,  TECH  INFO  FAC 
P.O.  BOX  8757,  ATTN:  ACQ  BR  I 


BALTIMORE/WASHINGTON  INTL  AIRPORT 
MARYLAND  21240 


COMMANDER 
DEFENSE  DOCU  CEN 

ATTN:  DDC-TCA  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 

METALS  $ CERAMICS  INFO  CEN 
BATTELLE  COLUMBUS  LAB  1 

505  KING  AVE 
COLUMBUS,  OHIO  43201 

MPDC  1 

13919  W.  BAY  SHORE  DR. 

TRAVERSE  CITY,  MI  49684 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  BENET  WEAPONS  LABORATORY, 

DRDAR- LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


